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Abstract 

We discuss the concept of natural Poisson bivectors, which allows us to consider the over- 
whelming majority of known integrable systems on the sphere in framework of bi-Hamiltonian 
geometry. 



1 Introduction 

The Hamilton- Jacobi theory seems to be one of the most powerful methods of investigation the dy- 
namics of mechanical (holonomic and nonholonomic) and control systems. Besides its fundamental 
aspects such as its relation to the action integral and generating functions of symplectic maps, the 
theory is known to be very useful in integrating the Hamilton equations using the variables separa- 
tion technique. The milestones of this technique include the works of Stackel, Levi-Civita, Eisenhart, 
Woodhouse, Kalnins, Miller, Benenti and others. The majority of results was obtained for a very 
special class of integrable systems, important from the physical point of view, namely for the systems 
with quadratic in momenta integrals of motion. The Kowalevski, Chaplygin and Goryachev results on 
separation of variables for the systems with higher order integrals of motion missed out of this scheme. 

Bi-Hamiltonian structures can be seen as a dual formulation of integrability and separability, in 
the sense that they substitute a hierarchy of compatible Poisson structures to the hierarchy of functions 
in involution, which may be treated either as integrals of motion or as variables of separation for some 
dynamical system. The Eisenhart-Benenti theory was embedded into the bi-Hamiltonian set-up using 
the lifting of the conformal Killing tensor that lies at the heart of Benenti's construction [U [15]. The 
concept of natural Poisson bivectors allows us to generalize this construction and to study systems 
with quadratic and higher order integrals of motion in framework of a single theory [311. 

The aim of this note is to bring together all the known examples of natural Poisson bivectors on 
the sphere, because a good example is the best sermon. Some of these Poisson bivectors have been 
obtained and presented earlier in different coordinate systems and notations. Here we propose the 
unified description of this known and few new bivectors using so-called geodesic H and potential A 
matrices |31j . In some sense we propose new form for the old content and believe that this unification 
is a first step to the geometric analysis of various natural systems on the sphere, which reveals what 
they have in common and indicates the most suitable strategy to obtain and to analyze their solutions. 

The corresponding integrable natural systems on two-dimensional unit sphere §^ are related to 
rigid body dynamics. In order to describe these systems we will use the angular momentum vector 
J = (Ji, J2, J3) and the Poisson vector x = {xi,X2,x^) in a moving frame of coordinates attached to 
the principal axes of inertia [1] . The Poisson brackets between these variables 

may be associated to the Lie-Poisson algebra of the three-dimensional Euclidean algebra e(3) with two 
Casimir elements 

3 3 
Ci = = C2 = (a;,J)=^XfeJfe. (1.2) 

fc=l k=l 

Below we always put C2 — 0. 
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As usual all the results are presented up to the linear canonical transformations, which consist of 
rotations 

X ^ aU X , J U J , 
where a is an arbitrary parameter and U is an orthogonal constant matrix, and shifts 

X X , J ^ J + S X , 

where S is an arbitrary 3x3 skew-symmetric constant matrix [?J [IB] . 

If the square integral of motion C2 — {x, J) is equal to zero, rigid body dynamics may be restricted 
on the unit sphere §^ and we can use standard spherical coordinate system on it's cotangent bundle 

xi — siii(j)sm9, X2 = cos(/)sin0, 0:3 = cos ^, 

(1.3) 

sin cos 6* cos cos 6^ 

Ji = ^-7 — P4, - coscppe , J2 = — — p^+smcppe, J3 ^ -Pci> ■ 

sm u sm u 

We use these variables in order to determine and classify the natural Poisson bivectors on T*SP up to 
the point canonical transformations. 

As far as the organization of this paper is concerned, in Section 2 we briefly introduce the notions 
of bi-Hamiltonian geometry relevant for subsequent sections. In particular, we discuss the concept of 
natural Poisson bivectors on cotangent bundles to Riemannian manifolds, which allows us to generalize 
classical Eisenhart-Benenti theory. In Section 3 we discuss the bi-Hamiltonian classification of bi- 
integrable systems on the sphere. Section 4 is devoted to the separable natural systems coming from 
auxiliary bi-Hamiltonian systems. 



2 Some issues in the geometry of bi-Hamiltonian manifolds 

A bi-Hamiltonian manifold AI is a smooth manifold endowed with a pair of compatible Poisson bivectors 
P and P' such that 

[P,P']=0, [P'.P']^0, (2.1) 

where [., .] is the Schouten bracket. This means that every linear combination of P and P' is still a 
Poisson bivector. 

If P is invertible Poisson bivector on M, one can introduce the so-called Nijenhuis operator (or 
hereditary, or recursion) 

N = P'P-^ . (2.2) 

If N has, at every point, the maximal number of different functionally independent eigenvalues 
Ml, . . . , u„, then M is said to be a regular bi-Hamiltonian manifold. 



2.1 Bi-integrable systems 

Let us consider a family of bi-integrable systems for which there are functionally independent integrals 
of motion Hi , . . . , Hn in the bi-involution 

{H,,H^}^{H,,Hjy ^0, z,j = l,...,n, (2.3) 

with respect to a pair of compatible Poisson brackets {., .} and {., .}' defined by P and P' . There are 
three known distinct constructions of bi-integrable systems, see |3l^ . 

Firstly, if M is a regular bi-Hamiltonian manifold endowed with invertible Poisson bivector P, 
then we can construct recursion operator N (j2.2p and, as usual, functions 

Hfc^^triV^- (2.4) 
form a bi-Hamiltonian hierarchy on M, i.e. the Lenard relations hold 

P'd-Hk = PdUk+i , for all A: > 1 . 
Using these relations we can get all the integrals of motion starting with the Hamilton function Hi. 



2 



Remark 1 The natural obstacle for existence of the bi-Hamiltonian systems is discussed in [5]. For- 
tunately, we can use these rare bi-Hamiltonian systems (natural or non-natural) as auxiliary systems 
for the construction of an infinite family of non bi-Hamiltonian separable systems. 

Namely, second special but more fundamental construction of integrable systems was originally 
formulated by Jacobi when he invented elliptic coordinates and successfully applied them to solve 
several important mechanical problems: "The main difficulty in integrating a given differential equation 
lies in introducing convenient variables, which there is no rule for finding. Therefore, we must travel 
the reverse path and after finding some notable substitution, look for problems to which it can be 
successfully applied". 

In framework of the Jacobi method we consider Hi (|2.4[) as constants of motion for an auxiliary 
bi-Hamiltonian system on the regular bi-Hamiltonian manifold M and treat functionally independent 
eigenvalues uj of N 

B{X) = ( det(iV - AI) ) = (A - ui)(A - U2) • • • (A - w„) , (2.5) 

as " convenient variables" for an ifinite family of separable bi- integrable systems associated with various 
separated relations 



^iiui,Pu,,Hi, . . . , Hn) ^ , i = with det 



^ . (2.6) 



Here u = (wi, . . . , w„) and Pu — {pm , ■ • ■ ,Pu„) are canonical variables of separation 

{ui,puj} = Sij and = S^j . (2.7) 

The Poisson brackets (j2.7|) entail that solutions Hi, ... , of the separated relations (|2.6p are func- 
tionally independent integrals of motion in the bi-involution (j2.3|) . see |26| . Of course, this construction 
will be justified only if we are capable to obtain separable Hamilton functions H — Hi, which have 
natural form in initial {p, q) variables (j2.10p . 

The third construction of integrals of motion in bi-involution on irregular bi-Hamiltonian manifolds 
is discussed in [HI [31] . In this case polynomial integrals of motion H2 , . ■ . , Hn are solutions of the 
following equations for the given Hamiltonian Hi 

P'dHi = XkPdlnHk, k> 1, Xk € R , 



which replace the usual Lenard relations (j2.4p . If this equations have many different functionally 
independent solutions labeled by different Xk , then we obtain so-called superintegrable systems [El [31] . 

2.2 Bi-Hamiltonian structures on cotangent bundles 

According to |32j a torsionless (1,1) tensor field L on a smooth manifold Q gives rise to a (second) 
Poisson structure on the cotangent space M = T*Q, compatible with the canonical one. 

Let 9 be the Liouville 1-form on T*Q and uj = dO the standard symplectic 2-form on T*Q, whose 
associated Poisson bivector will be denoted with P. If we choose some local coordinates q = {qi, . . . , qn) 
on Q and the corresponding symplectic coordinates {q,p) — {qi, . . . ,qn,Pi, ■ . . ,Pn) on T*Q then we 
get the following local expressions 

9 = pidqi + . . .pndqn , and ^ ^ ^ \ Q ) • (^"^^ 
Using a torsionless tensor field L one can deform 6* to a 1-form 9' and P to bivector P': 

( \ 

n 

9' = LijPidqj, and P' = 



V " ^1 V dqj dq, J J 

The vanishing of L torsion entails that P' (|2.9p is a Poisson bivector compatible with P. 
Let us consider natural integrable by Liouville system on Q. 



(2.9) 
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Definition 1 The natural Hamilton function 



1,3 = 1 



(2.10) 



is the sum of the geodesic Hamiltonian T defined by metric tensor g{qi, . . . , g„) and potential energy 
V{qi, . . . ,(7„) on Q. 

If the corresponding Hamilton- Jacobi equation is separable in orthogonal coordinate system {u,pu) 
on configurational space Q, then in framework of the Eisenhart-Benenti theory the second Poisson 
bivector P' (|2.9p is defined by a conformal Killing tensor L of gradient type on Q with pointwise 
simple eigenvalues associated with the metric g{qi, . . . , g„), see [TJ [51 [31 151IT3]. 

According to Kowalevski [T7j and Chaplygin [7], separation of variables for integrable systems 
with higher order integrals of motion involves generic canonical transformation of the whole phase 
space. Definition (|2.10p of the natural Hamiltonian and metric tensor g(gi, . . . ,qn) is non-invariant 
with respect to arbitrary canonical transformations of coordinates on T*Q 

qi q[ = Mp, q) , Pi Pi = 9i{p, q) ■ 

In the situation, when habitual objects (geodesic, metric, potential) lose their geometric sense and 
remaining invariant equation (j2.ip has apriority infinite many solutions, notion of the natural Poisson 
bivectors on T*Q became de-facto very useful practical tool for the calculation of variables of separation 

[HIIHllMllMliSIlIM]- 

Definition 2 The natural Poisson bivector P' on T*Q is a sum of the geodesic Poisson bivector P^ 
compatible with P 

[P,P^] = [P^,P^]-0, 



and the potential part defined by a torsionless (1,1) tensor field A(qi, ... ,qn) on Q 



P' = Pfr 



( 







k=l 



\ 



dqj dqi 



Pk 



(2.11) 



(2.12) 



In fact, here we simple assume that bi-integrability of the geodesic motion is a necessary condition for 
bi-integrability in generic case at V 0. 

Throughout this paper geodesic bivector PJp is defined by n x n matrix Tl{qi, . . . , qn,Pi, ■ ■ ■ ,Pn) 
and functions x, y and z on T*Q 



P' 



k=l 



V 



dpi 



-H, 



dpj 



Hi 



\ 



dn 



k] 



k=l 



^kip) 



(2.13) 



up to the point transformations. In this case the corresponding Poisson bracket {., .}' looks like 



{qz^Pj}' = n„ +A,j , {qz,qj}' = E ~ Yikiq)^^'''' 

k=l 



k=l 



OA 



ki 



dA 



kj 



dq, dq^ 



Pk 



E 

k=l 



dp, 



\ dq. 



dqi 



dpj 
zk{p)- 



In fact, functions x, y and z are completely determined by the matrix H via compatibility conditions 
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We can add various integrable potentials V to the given geodesic Hamiltonian T in order to get 
integrable natural Hamiltonians (j2.10p . In similar manner we can add different compatible potential 
matrices A to the given geodesic matrix 11 in order to get natural Poisson bivectors P' (|2.12p compatible 
with the canonical bivector P. 

Remark 2 We have to underline that this definition of natural Poisson bivectors is the useful anzats 
rather than rigorous mathematical definition. It is an obvious sequence of non-invariant definition 
of the natural Hamiltonian with respect to transformations of the whole phase space. We hope that 
further inquiry of geometric relations between n x n metric matrix g, potential matrix A and geodesic 
matrix II on T* Q allows us to get more invariant and rigorous definition of these objects. 

Remark 3 In term of variables of separation 11 = and A — diag(iti, . . . so we have usual 
invariant construction of Turiel [32] . The main problem is how to rewrite this invariant theory in term 
of initial physical variables. 

Remark 4 We suppose that (j2.13|) is a special form of P^. Other form of P^ on the generic symplectic 
leaves of e*(3) for the Steklov-Lyapunov system at C2 7^ will be presented in the forthcoming 
publication. 

3 Special natural Poisson bivectors on the sphere 

The standard Laplace method for the direct search of integrable systems may be applied to the search 
of the natural bivectors P' too. 

Firstly we stint ourselves by a family of natural Poisson bivectors (|2.12[) with geodesic part ()2.13|) . 
Then, it is easy to see that the geodesic Hamiltonian 

n 

on the cotangent bundle T*Q is the second order homogeneous polynomial in momenta, so we assume 
that entries of H are the similar homogeneous polynomials 

n 

E (3-1) 

k,m—l 

up to canonical transformations pk Pk + fk{<lk)- On two-dimensional unit sphere Q = we use 
spherical coordinates (|1.3|) such that 

a = (91,92) = ((^^',^') and p^{pi,p2)^{p4„pe)- (3.2) 

At the third step we introduce a family of partial solutions for which all the entries of (|2.13p are 
independent on variable 0, i.e. at 

c-r(9) = 4'"(^)' x,fe(0,0)=x,fc(0), y,k{4>,e)=y,k{e). (3.3) 
It looks like reasonable assumption because the geodesic Hamiltonian T 

T = ai Ji -I- a2 J| -I- as Jg = ^03 - cot^ 6'(ai sin^ + a2 cos^ 0)^p| , 

(3.4) 

— sin2(/)Cot^?(ai — a2)pBP4> + (oi cos^ + 02 sin^ (/))pg 

is independent on variable at ai — a2. If are constants it means that two diagonal elements of 
inertia tensor of the body a^^ =0-2^ are equal to each other and we discuss symmetric rigid body [3]. 

Due to the special form of (12.13^ and additional assumptions (|3. 1113. 31) . equations (12. lip de- 
compose on the subsystem of equations for c^j^{q), subsystem of equations for Zk{p),c'^j^{q) and third 
subsystem of equations for Xfc(g), yfc(g), c^™(q), which can be partially solved independently to each 
other. 
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Proposition 1 // assumptions !l3.1Vi3.3\) hold, then subsystem of equations for the functions Zk{p) 
coming in h2.11]) has three families of solutions 



Case 1. Ylij = 0; 

Case 2. Zl = , Z2 = ; (3 5) 

Case 3. zi = — , Z2 — — . 

3 ' 3 

This proposition gives only ttie necessary conditions. Of course, there remain complementary equations 
on the other functions c*™(0), Xjk{0) and yjk{d) which have to be solved in the sequel. 

At the first case = and we can immediately look for compatible potential part A((/), 9) and 
the variables of separation M1.2 (|2.5I) . which are related with initial variables by the point canonical 
transformations 

Ui ^ fii(t>,0), ^ gi{4>,0)p^ + hi{(j),e)pe . (3.6) 

As a consequence, the geodesic Hamiltonian is a second order homogeneous polynomial in physical and 
separated momenta and the theory of projectively equivalent metrics in classical differential geometry 
study essentially the same object [3]. 

Proposition 2 In second case generic solution of is parameterized by six functions g, h and 

one parameter 7 = 0, 1; 

/ Wl 91 {0)pI + 92 {0)P4. Pe + 53 {0)pI \ 

n= (3.7) 
V hi {0)pI + h2 {e)p^ pe + /13 {e)pl ) 

up to the point transformations pk — > akPi + PkP2 ■ 

As above it is only necessary condition and functions 5, h from p.7p . together with functions x, y from 
()2.13p . are solutions of the remaining six non-linear differential equations in p. lip . 

Proposition 3 In third case generic solution of h2.11\) is parameterized by nine functions f,g,h and 
one parameter 7 = 0, 1.' 



/ fiiO)pl + f2{0)p^Pe + fmvl gmPl + 92iO)P4,Pe + 93{0)pI 



n 



"0 

3/2 



V ^/2(%'^ + 2/3(^^)p0Pe+7(/3(^) + /i3(e)) pI h,{e)pl + hMP<pPe + h^{0)p' 
up to the point transformations pk — > OLkPi + PkP2 ■ 



(3.8) 



,2 



Functions /, 5, h from (|3.8p . together with functions x, y from (|2.13p . are solutions of the remaining 19 
non-linear differential equations in (|2.1ip . 

Matrices p.7p and p.8p were obtained as solutions of the subsystem of algebraic and linear differ- 
ential equations for cfp{9), which has an unambiguous solution. The remaining functions satisfy to the 
complementary overdetermined subsystem of nonlinear PDE's, which have many distinct particular 
solutions. 

In both cases (|3.7I) and (13.81) we can get a complete classification of these particular solutions and 
of the corresponding bi-Hamiltonian systems (|2.4I) . Classification of separable bi-integrable systems 
demands additional assumptions on the form of the separated relations. 

3.1 Case 2 - classification of natural bi-Hamiltonian systems 

Let us briefly discuss a procedure of classification of the natural bi-Hamiltonian systems associated 
with natural Poisson bivector (12. 12112.151) defined by the geodesic matrix 11 (|3.7I) . 
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If h2{9) = in p.7|) . then six differential equations coming in (|2.1ip have four distinct solutions; 
among them we pick out solution defined by the following matrix 



n 



h',{9)F 



F2 



F = tan a 



d9 



If fli = 02 = const, then we can put /i3(^?) = 7 = 1 without loss of generality and obtain 



n = 



yi2{e) 



2ayi22{0) — cosQ!0sinQ!0 



{I + t&T? aO) pI + p\ 
The corresponding geodesic Hamiltonian (|2.4p is equal to 



At a = 1 matrix 11 
/ 

A = 



(3.9) 



r = i tr iV = tr n = (2 + tan^ ae)pl + pi . 
is consistent only with the following potential matrix 
/(0) 9{^,0) 



where 



V 2cos6l 



/(0) 



.9(0,^) 



a cot (p 



2cos2(l){2cos^9 + l)g{(f),d) f{(j)) 



sin 20 sin 29 



SIT? 6 cos^ I 



cos^ 9 



a tan^ 9 



(3.10) 



2dcos2(?!)(2cos2(/)- 3) 



sm 



9(^,9) 



2d ail? 9 sin 2(j) 
cos 9 



So, bi-Hamiltonian system associated with 11 p.9p and A p.lOp has the following Hamilton function 
(EH) 



-Hi 



= T 



a{{xl + xl) - xlixj - xD) {l + xl){xl+xl){bxl + c{xl+xl)) 



2 2 



2 2 2 

rrt-^ /y-^ 



2d{xl + xl + 2xlxl){{xl + xl) - xlixj - x^)) 
{x\ + x'^)x\x'i^ 



Second integral of motion I-L2 (|2.4p is a fourth order polynomial in momenta. This integrable system, 
to the best of our knowledge, has not been considered in literature yet. 

In similar manner we can get a complete classification of natural bi-Hamiltonian systems associated 
with matrices (I3.7P and (I3.8p . 

3.2 Case 3 - one possible generalization 

Non-invariant assumptions (13.113. 3p depend on a choice of coordinate system and we miss a lot of 
another solutions of (|2.ip . which may be interesting in applications. 

One of the possible generalizations consists in the application of multiplicative separable functions 
in (O) 



a'ypi^)b1pi9), 



and similar for x, y. For instance, geodesic matrix 
(sin 9pg + i cos 9 p^ ) 





2 a^p0(sin6'p6i + icos^p^) 



H 



sm 




a e 



(3.11) 
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gives rise to the natural Poisson bivector P' at 

i P<t> Pe 

It is easy to prove that integrals of motion for the Lagrange top (|4.2p are in involution with respect to 
the corresponding Poisson bracket {., .}'. 

Remark 5 Bivector P!^ (|2.13p associated with 11 (|3.11l) has a natural counterpart on the generic 
symplectic leaves of the Lie algebra e*(3) at (a;, J) ^ 0. 



3.3 Case 3 - three-dimensional sphere 

On the three and four dimensional spheres endowed with the standard spherical coordinates there are 
the same three families of solutions p.Sp . It means that factor 1/3 in (13.51) is independent on dimension 
of the sphere. 

For instance, \i q — (<p, ip, , 6) and p ~ {p^,p^,pg) are the standard spherical coordinates on T*S^, 

then at zi. = — matrices 
3 



Hi 



pI 2p^p^ 
Pl + fpl + j^ 
■ Ptppe 



Ho 



Pi 






f,2 

^p^tp^ 
F 

-276""''' 



2ff"\ \ 

- -pr ) P'pP'^ 

f 2 CtP 2 



/ = /(^): 



(3.12) 



Pi^pe 



2(6"''' + l)i 
2l3e-''^{e''^ + l)^p^P9 
F - 47(6""''' + 1)pI 



where F = (e"''' + l)p^ + Pe'^'^ie"'^ + if p^^ + 76""''' p^, determine geodesic Poisson bivectors (^1^ 
and geodesic Hamiltonians (|2.4p 



Ti = 3p^ + Y p.^ + -jjY- Pe : 



r2 



(26"'> + 3)p^ + ^ 



■4-27(2+ — 



Then we can calculate compatible potential matrices A1.2 depending on coordinates (0,-0,^) and the 
corresponding integrable potentials Vi^2- The corresponding integrals of motion 7^2,3 (12.41) are the 
fourth and sixth order polynomials in momenta, respectively. 

So, using notion of the natural Poisson bivectors we can produce a lot of abstract mathematical 
examples of bi-Hamiltonian system on the sphere. The main problems are how to select physically 
interesting bi-Hamiltonian systems and how to construct significant separable systems from the non- 
physical auxiliary bi-Hamiltonian systems. 



4 Separable bi-integrable systems 

In this Section we present matrices H and A for the following well-known separable systems on the 
sphere 

• Case 1 - Lagrange top, Neumann system and systems separable in the elliptic coordinates; 

• Case 2 - Goryachev system, Matvcev-Dullin system, Kowalevsky top, Chaplygin system; 

• Case 3 - Goryachev-Chaplygin top, Sokolov system, Kowalevsky- Goryachev- Chaplygin gyrostat; 
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which may be natively embedded into the proposed scheme as separable bi-integrable systems. Some 
new mathematical generalizations of these systems and new separation of known systems are collateral 
results for this activity. 

In framework of the Jacobi methods one gets integrals of motion Hi, ... , Hn as solutions of the 
separated relations (|2.6p . Of course, variables of separation and separated relations could have the 
singular points. So, the standard problem is the rigorous determination of domain where variables of 
separation and integrals of motion are well defined, see the Jacobi definition of the elliptic coordinates. 

Our main purpose is to discuss natural Poisson bivectors and, therefore, we do not comment this 
huge and comphcated part of the work here, see for example [H [71 IH [TOl [171 [35] and references within. 



4.1 Case 1 - Lagrange top 

If the spherical coordinates (j), 9 (I1.3P are variables of separation, one gets the simplest natural Poisson 
bivector P' ([^1^ at 

n = 0, and A= ^ ^ ^ . (4.1) 
The auxiliary bi-Hamiltonian system is trivial 

On the other hand, substituting variables of separation ui ~ (j) and U2 ^ into the separated relations 
/ cos^ \ 

$1= [a+ — ^ i?2-i?i+]3e + focos6i = 0, *2 = - = , 
V sin 6* / 

one gets integrals of motion for the Lagrange top in rotating frame 

Hi = Jf + J'i + aJl + bx3, H2 = Ji, a, 5 e M , (4.2) 

More complicated natural bi-vector P' obtained from matrix gives rise to another variables of 

separation for this system. 

Remark 6 According to "FP^ bivector P' (|2.12p associated with A (j4.1l) admits extension from cotan- 
gent bundle T*§^ to the symplectic leaves of the Lie algebra e*(3) at {x, J) ^ 0. 



4.2 Case 1 - Neumann system 

Let us put = in ()2.12p and consider some particular solution P' of the equations (|2.ip defined by 
the following non-symmetric matrix 



A 



V 



ai cos -f a2 sm 
(fli — 02) sin 



(fli — 02) sin 20 cos Q 



sin0 



cos d sin Q az sin B + (ai sin (/) + 02 cos (/)) cos d 



(4.3) 



with three arbitrary parameters G R. As above, the auxiliary bi-Hamiltonian system has trivial 
integrals of motion 'H.k (12. 4p . which are functions only on the configurational space S^. 

On the other hand, coordinates of separation Uj (|2.5p are the standard elliptic coordinates on the 
sphere 



(A - ui)(A - U2) 



A - ai A -02 A -03 (A - ai)(A - a2)(A - as) 
By substituting these variables in the separated relations 

u^Hl - H2- 4(ai - Ui){a2 - Ui)(a3 - Mj) vt. + Ui{ui) = 0, 



(4.4) 



1,2, 
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one gets bi-integrable systems with quadratic in momenta integrals of motion 

Hi = Jl + Jl + J| + V{x) , H2 = aiJl + 02^2 + «3 J| + W{x) , 

which are in the bi-involution (|2.3p with respect to both Poisson brackets. Here V{x) and W{x) are 
easy calculated from the potentials J7i,2- For instance, if 

U{u) — u{u — fli — 02 — 03) , 

then one gets the Neumann system with the following integrals of motion 

Hi ~ Jf + J2 + Ji + aixi + 02X2 + 03X3 , 

(4.5) 

H2 = fll Jl + 02^1 + 03^3 - 0203X1 - 0103X2 - 0102X3 . 

Remark 7 Bivector P' (|2.12l) associated with A (|4.3p also satisfies equations (|2.ip at (x, J) ^ 0, but 
in this case we lose bi-involutivity (|2.3p of integrals of motion iJi,2 (|4.5p for the Clebsch system on the 
whole phase space e*(3). Of course, the corresponding elliptic coordinates on e*(3) remain variables 
of separation, but we can not get interesting natural Hamiltonians using these variables |25j . 

4.3 Case 2 - systems with cubic integral of motion 

At 7 = in p. 71) we have particular solution of the equations (|2.1ip defined by geodesic matrix 
2 {d9 ^ g{e) ) 



n = 



F = 



[9{9)pe ~ ih{0)p^) , i = V^, (4.6) 



V F J 

depending on arbitrary functions g{0) and h{0) and by functions 

n n 

This matrix 11 is consistent with the diagonal potential matrix 



A = a exp ( 10 - / ^ ( ° I . (4.7) 



1 



The corresponding bi-Hamiltonian systems (|2.4p are non-physical T = F and, therefore, we immedi- 
ately proceed to consideration of the coordinates of separation wi.2 = ^/W^ following to |34) . If we 
introduce polynomial 

6(A) = (A - i;i)(A - 1^2) = A^ - iVFA + Aij . 

instead of characteristic polynomial B{X) = (A — mi)(A — U2) = (A — vf){\ — ul) (|2.5p of recursion 
operator N, then it is easy to prove that 

fi-X \ X n J 

\d6 ip0 



where 



It entails that 



J 9{S) A 
p^. = A{X = Vj) , j = 1,2, 
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are canonically conjugated to Vj momenta and that the corresponding Poisson brackets read as 

{Vi,Pv^} = <5y , {vt,Pvjy = Sijvf . 

Now we have to substitute this family of variables of separation into the separated relations and try 
to get natural Hamiltonians. For instance, let us take 

g{0) = sm0f{e), h{e)=cosef{e), (4.8) 

substitute 

2i 

X = Vj ^j.= —VjPj, j = l,2, 

into the equation 

2 

$(A, fi) = fiHi + - - - 6A + ^ = 0, (4.9) 

A 

and solve a pair of the resulting equations with respect to Hi 2- If ai — 02 in the geodesic Hamiltonian 
p.4[) . then in this solution we have to put 

sm 

and we obtain integrals of motion for the Gorychev system on the sphere [12] 



Q -5 ' 2/3 2/3 

•Juri djr) 



H 2 J3 / 2 ^ 2 _L 8 2 b \ 1/3 4ia 



(4.10) 



For other separable natural bi-integrable systems from |23[ 134) we present Hamiltonians and functions 
g and h only. So, for the Goryachev-Chaplygin top [TTl [7] we have 

7J, = J,2 + j|+4j| + axi + A, g(^e)^^— /.(^)^ 37'^.-^ 

cos sm cos"^ 6 sm 



For the DuUin-Matveev system [9] with Hamiltonian 



Hi = j2 + j| + ( 1 + - ;/ '"' ) 



X% — \X\ \ 2 OL^l 



x^ + c 4(.X3 + c) V (2:3 + c) 1/2 a;3 + c 

geodesic matrix H ()4.6p and potential matrix A (|4.7p are defined by functions 

1 1 -2ccos6i-3cos2 6' 

^ ^ sm6l ^ ^ 2sin2 0(cos6l + c) 

For the system with the Hamiltonian 



12 2(x3 + |x|); ^ (a;3 + N|)5/6 (3:3 + |:e|)i/3 
bi-Hamiltonian structure is defined by functions 

_ (cOsg + l)V3 (cOSg + 1)^/^ 

^h^^ ' ^(^)--2(cos0)-l) ■ 

For the last system from [53] we have 

t2 , t2 , A3 3X3 \ j2 , , b 



16 8(X3 + |X|); ^ (a;3 + |x|)3/4 (x3 + |x|)l/2 
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and 

(COS0+ 1)1/2 (3cos0+l)(cos0+ 1)1/2 

sm 4 sm Q 

If P' is the linear in momenta Poisson bivcctor from [34^ . then our natural Poisson bivector is equal 
to P' = P'P-^P' . 

Remark 8 According to [33] , the Coryachev-Chaplygin, Chaplygin and DuUin-Matveev systems can 
be embedded into a family of integrable systems with cubic integral of motion. We suppose that bi- 
Hamiltonian structures for the Valent systems may be described by a suitable choice of the functions 
g{Q) and h{B) in dM]) and 

Remark 9 Another possible generalization consists of multiplication of matrix 11 (|4.6p on the functions 
depending on similar to p. lip . 

4.4 Case 2 - Kowalevski top and Chaplygin system 

Let us consider a geodesic bivector Plj, (|2.13p determined by the matrix 11 

/o ^ \ 

^=-^a —a " , aeK, (4.11) 

sm Q cos"' Q \ n 

V co^^Qv\ + ^^^Qvl I 

and by functions 

yi2(^) cos (?^sin Q + ax22(0) cos , zi^2 = . 
There is only one potential matrix consistent with 11 (|4.1ip 



/ a cos — b sin a(j) (a sin acj) — b cos a(j) \ cot \ 



A 



a,beR. (4.12) 



(4.13) 



y sin a(f> — b cos a(l)j ta.n 9 —a cos acj) + b sin acp j 

The corresponding coordinates of separation ui^2 (|2.5p are the roots of the polynomial 

2 pI sin^ +p'^ cos^ 9 (a cos a(j) ~ b sin acj)) (pi sin^ & + cos^ 9) 

B(\) = A — t: A — 

^ ' sin"6'cos2 

2 sin 9 (a sin acj) + b cos acj>)p^ pg 2 ,2 
sin"6'cos2 6l " ^ ' 

Following to [29l [30] we can introduce auxiliary polynomial 

sin9pg a sin acj> + b cos acj) sin 9( a cos acj} — b sin acj>) 

^ A) = A + p^ pe , 

a cos 9 a a cos 9 

such as 

{BiX),A{p)} - ^ ((/.2 _ ^2 - 62)B(A) - (A2 - - 62)B(^)) , {A{X), A{f,)} = . 
It entails that ^ 

are the canonically conjugated momenta satisfying to the Poisson brackets (|2.7p . At a = 2 these 
variables have been considered by Chaplygin [S]. 
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By substituting these variables of separation into a pair of the separated relations 

$1 - [ul -a"- b^)pl^ +Hi-H2=0, $2 - (^^2 - - 62)^2^ + + H2 = , 
one gets separable bi-integrable system with the Hamilton function 

Hi = p^- tan^ Opl + 2{acosa(t> + b cos acj)) cos" 6 , a eR. 
According to [3D] , at a = 1 using separated relations 



(4.14) 



one gets Hamilton function of the generalized Kowalevski top [T7] 



jjkotv = 2ffi = ( 1 - 



c+1 



(Jf + J|) + 2J| + 2aa;2 + 2bxi 
At a = 2 we can use another separated relations 



+ du = (4.15) 



(4.16) 



$(u,p„) = ((m^ _ a2 _ ^2^,^2 + + Fi - iJ2) ((u^ - - b'^)pl + cu + Hi + H2) + du ^ (4.17) 
in order to get Hamiltonian of the generalized Chaplygin system [5J [T3] 



H 



ch 



8H1 



1 



4c+ 1 



{Ji + J2) + 2 J3 - 2a(x? - x^) - 2bxiX2 



2d 



1 + 4c - 



(4.18) 



At c = —a 2 we have geodesic Hamiltonian T = Ji + J2 + 2 J| with the constant inertia tensor. 

Remark 10 By substituting these variables of separation into another separation relations we can 
obtain various mathematical generalizations of bi-integrable Hamiltonians (14.1414. 16I4.18P . 

4.5 Case 2 - spherical top and Chaplygin system 

At 7 = in p.7p we have a particular solution of the equations (j2.1ip defined by matrix 

a - sir? 6 



( 



H 



cos^ 9 sin^ 



:P<kPe 



\ 



a e 



(4.19) 



sm 



and functions 



. ^ ^ 2acos2 6l 

yi2 = sm 6* cos 6* + , „ X22 

sm U — a 

In this case coordinates of separation ui^2 (|2.5p are equal to 



Zfc = . 



Ul=P4,: 

SO that conjugated momenta read as 



U2 



apl (sin^ 9 - a)pl 



sm^9 



arctan 



t tan 9 



sm9 cos 9 arctan 



sm Upg 



I a cos^ 07)2 — sin^ ^(sin^ i 



2P0 



P«2 



2^ /a cos 



2 fl„2 



sm W sm 



a)pl 
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By substituting these variables of separation into the separated relations 

$1 = - = , $2 = ai^i - ^^2 (l - (a - 1) tan2(2p„, ) + af{e) = 
where 



arccos * 



'u2-aH2 a(if| - U2)(l - cos 4p„2^/uJ) 



y U2 2U2 

one gets generalized Lagrange top with integrals of motion 



Hi = Jl + J| + 4 + f{x3) , H2 = J3 



(4.20) 



Ml sin (4p„j 



H2- Hi+ui=0, 



Other separated relations 



$2(^1, Pui) affi - 1*2(1 - (a - 1) tan^(2p„2 0IJ)^ 

give rise to integrals of motion for the spherical top 

Hi = T = Jl + J2 + , H2 = J1J2J3 ■ 

There are only two potential matrices compatible with 11 (|4.19p 



(4.21) 



(4.22) 



\ 



\ 2 sin 6 cos 6 



and 



A(2 



sm f 



v 



a sin 20 + 6 cos 20 



- {a — sin^ 9){a cos 20—6 sin 20) 



a sm f 



- {a — sin^ &) (a cos 20 — 6 sin 20) 



(a - 2 sin^ 0) (a sin 20 + 6 cos 20) 



a cos t/ a. 
In the first case the auxiliary bi-Hamiltonian system with the Hamilton function (|2.4p 

a- r 



H*^) = ( 1 



Jt + J2^ + 2Ji + / 



/ Xl_ 

\X2 



1 



(4.23) 



is a deformation of the geodesic Hamiltonian for the Kowalevski top at a = 1 and / = 0. By 
substituting the corresponding coordinates of separation 



wi='"i + /(0), j5„i=p„i-i 



dx 



^ a/(0) „ 

U2 = "2 + . 2/1 ' -P"2 = -P^i2 ) 

sm 6* 



into $1 = Ml — i/2 = and the second separated relation $2 in (|4.2ip . one gets a generalization of the 
spherical top defined by the following integrals of motion 



/ 



Hi = J]^ + J2 + 0/3 



H-i — J-x 



X2 



37-1 I 3Jn I 3jQ 



In the second case matrices 11 (j4.19p and A^^^ give rise to the auxiliary bi-Hamiltonian system with 
the Hamilton function 



= ( 1 



a- 1 



^ (ji + jt) + 2J| + 4a"ia2;ia;2 - 2a^H{x\ - a;^) 



(4.24) 



It is a new deformation of the well-known Chaplygin system [B] . 
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Remark 11 According to [21], there is a non-canonical map, which relates integrals of motion (|4.22p 
with integrals of motion for the Gaffet system [TU] 

TIT _ x2 , t2 , t2 1 H- - 7 r 7 , X2X3.J1+X1X3.J2+X1X2J3. 

(XiX2X3)'^'-^ [XiX2X3y'-^ 

In order to describe the bi-Hamiltonian structure for the Gaffet system we have to use additional 
non-point transformation of the standard spherical coordinates, which changes the form of (|2.13p 
in initial variables. This bi-Hamiltonian structure will be discussed in the forthcoming publication. 

4.6 Case 3 - Goryachev-Chaplygin top and Sokolov system 

At 7 = in (|3.8p equations (j2.1ip have a particular solution (|2.13p defined by the following 
symmetric matrix 

+ P0 (4 + 3 cot^ aO) 2p^ pg 

n = ( I , a e M, (4.25) 

2p0P9 pl-plcot^aO 



and by the functions 

cos aO sin a9 



Pk 



X22 = yi2 = , Zfc = o 

a 6 



There is only one potential matrix compatible with 11 ()4.25p 



A = 



1 



cos^ ad \ 1 

The corresponding auxiliary bi-Hamiltonian system is defined by the Hamilton function ()2.4p 

i^i = (2 + cot^ ae)pl +pI + — ^ . 
2 ^ cos^ aO 

If a = 1, we have a deformation of the geodesic Hamiltonian for the Kowalevski top [17] 

1 22 2 

-Hi — Ji + J2 + 2 J3 4 — 2 ■ 

This auxiliary bi-Hamiltonian system gives rise to the variables of separation ui^2 (12.51) 

2 



..2=|p,±\/-^+I.^ + ^j -[js±^JnJi + Ji + ^J , a = l. 

At a = these coordinates were found in [7 . By substituting the generalized Chaplygin variables 

1 / \ \n{vf 2 ~ o,) 
vi,2 = V"i.2 : = -^ln(^ui,2(ia;i - a;2) - (iJi - J2)a;3j H |t , (4.27) 

into the separated relations 

^i,2iv,pv) = Hiv + H2 + b\/ v'^ - asi-a2py - - cv"^ = , « = x'i,2, = Pi>i,2 > 
one gets integrals of motion for the generalized Goryachev-Chaplygin gyrostat [TTJ [^ 

Hi = Jl + ,4 + 4 J| + 2cJ3 + bxi + -^ 

X3 

H2 = (2 J3 + C) (jl + J2 + - &:E3 Jl . 
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By substituting the same variables (|4.27|) into the following separated relations 



$1,2(^1,2, Pt;i, 2) ^ Hi±H2 + b^vj ^ - a sin2p„i 2 - vl ^ - cvi^2 = 0, (4.28) 
we obtain the generalized Sokolov system [5D] defined by integrals of motion 
Hi = Jf + JI + 2JI + ch + h{J3Xi- xzJi) + , 



a 

3 ^ 32 ' 
■^3 



H2 = (2J3 +C + bxi^ ^1 J2 + J| + J, 

up to the canonical transformation discussed in [TB]. 

4.7 Case 3 - Kowalevski-Goryachev-Chaplygin gyrostat 

The geodesic matrix 11 (I4.25P for the Goryachev-Chaplygin top may be deformed 

^ ^ cos a9 



n = n + /3 



sin^ a9 ' 
V 



(4.29) 



if 

yii(6') = X2i(6l) - — , 2/12(6*) = r^— -X22(6l) : -, Zfc = — . 

2a sm a& a sm at) 3 

Remark 12 In the r-matrix formalism transition from matrix (|4.25p to the matrix (|4.29p generates 
transition from the quadratic Sklyanin bracket to the so-called reflection equation algebra [331 HH] • 

In generic case matrix 11 (|4.29l) is compatible with the potential matrix 

_w / cos^a0~A _£££i^ \ ^ . 2 ^ / 1 5 \ 

^ = 2 — + — ^r-^ (4-30) 

sin ae y 4 cos gfl singe ^^^2 J COS a9 ^ ^ ^ j 

The corresponding auxiliary bi-Hamiltonian system is defined by the Hamiltonian 

1 , 9 s 9 9 a(cos^ q;0 — 2)e^~^ hsiv? aO 
-Hi = 2 + cot^ a0)pl +pI + ^ r^-^ + — ^ • 

2 ^ sm aO cos"^ aO 

So, at a = 1 we have another deformation of the geodesic Hamiltonian for the Kowalevski top [17] 
l-iy _ r2 , t2 , 72 a{xl+xl + l) _^_4a,ot.n^.i/x2) ^ b{xl+xl) 



In this case description of the variables of separation and the corresponding bi-integrable system is an 
open problem. 

At /3 = ±2 i there is one more particular potential matrix compatible with 11 (|4.29p 

(± i sin aO cos aO \ 
(4.31) 
; 

In this particular case we can substitute the coordinates of separation ui^2 (I2.5P and the corresponding 
momenta 2 i'^to the separated relations defined by 



$('u,p„) = + Hiu^ + H2U^ + a+ ^/b{u)sm2pu = 0, (4.32) 
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and obtain integrals of motion for the Kowalevski-Goryachev-Chaplygin gyrostat with the following 
Hamilton function 



Hi = Jf + JI + 2 J| + 2ciJ3 + C2Xi + c^i^xl ~ xl) + 



C4 



(4.33) 



see |6l[T3l|T7l|33- Here b{u) (|4.32p is a special polynomial of eight order in u with coefficients depending 
on a and c^, see details in |22) . 

Remark 13 In this case in order to get the conjugated momenta Pm ^ and the separated relation we 
used the Lax matrices and the reflection equation algebra, that drastically simplified all the calculations. 

Remark 14 For the systems with quartic integral of motion from [53] the natural Poisson bivector 
may be obtained using deformation of the matrix (|4.29p similar to (|3.11l) . 

4.8 Case 2 - deformations of the Kowalevski top and Chaplygin systems 

Let us consider trivial canonical transformation 

VB^VB^!{Q). (4.34) 

which preserved canonical Poisson bivector P (|2.8I) . This mapping shifts the natural Poisson bivector 
P' (|2.12p associated with matrices H (|4.11l) and A (|4.12p by the rule 



/ 



P' 



m 







a cos^ 9 — 1 cot ( 



asm 



* 



* * 



cos W sm 



a 

,a-2 



P9 



sin" 9( a sin a(f> + b cos aa 



In' g 



■ (^sin 6 cos 9 In' g — 1 



where 



m 



2f{9) sin^ 



and 



cos^ 



In'g 



1 dg{9) 
9{0) de 



The Poisson bivector P' gives rise to the " shifted" variables of separation 



'lpe^P8+/(e) 



Pu - Pn\pg-^pg+f(e) ■ 



(4.35) 



If we substitute these variables of separation into the old separated relations (|4.15p and (|4.17p one gets 
non-natural Hamiltonians, which are related to the old Hamiltonians (|4.16p and (I4.18P by canonical 
transformation (|4.34p . 

In order to get new natural Hamiltonians we have to appropriately modify the separated relations. 
For instance, let us take 

V? tan"-i 9 
cos" 

At a = 1 by substituting variables of separation (|4.35l) into the new separated relations 

$ = $-/3i7i+/32 + y^(M2„a2-62)p„^ ?i = Mi,2, Pu^Pu^.,, (4.36) 
where $ is given by (|4.15p . one gets generalization of the Hamilton function (14.161) 



c+1 



d 13 



17 



At a = 2 the "shifted" separated relations 

$ = $ + V^(u^ -a^ - \?)i>u , u = ui,2 , Pu = , (4-37) 
where <& is given by (|4.17l) , yield similar generalization of the Hamiltonian (|4.18p 

i?- = (l - 1^-^) {Jl + Jl) + 2Ji 2ai4 - 4) ~ 2bx,., - + " 4) ■ 

These Hamiltonians at c = — and another Hamiltonians associated with various functions f{9) 
may be found in [55] , 

The separability of these systems, to the best of our knowledge, has not been considered in liter- 
ature yet. In tboth cases equations of motion are linearized on the two copies of the non-hyperelliptic 
curves of genus three defined by (|4.36[) and (|4.37p . We do not know how to solve the corresponding 
Abel-Jacobi equations as yet. 

Remark 15 Other natural Poisson bivectors studied in the previous Sections may be shifted on the 
similar linear in momenta terms. As above, it allows us to get various generalizations of the considered 
bi-integrable systems. 

5 Conclusion 

We proved that almost all known integrable systems on the two-dimensional unit sphere § may be 
studied in the framework of a single theory of natural Poisson bivectors. It is an experimental fact 
supported by all the know constructions of the variables of separation on the sphere. We try to draw 
attention to this experimental fact in order to find suitable geometric explanation of this phenomenon. 
So, this collection of examples may be helpful for investigations of the invariant geometric properties 
of metric g, geodesic 11 and potential A matrices as objects on the whole phase space, which allows 
us to obviate a necessity of the direct solutions of the equations ()2.1|2.1ip and p.3|) . Moreover, it can 
possibly be a suitable step towards the construction of Poisson bivectors on more generic symplectic 
and Poisson manifolds. 
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